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Abstract. Interaction of impulsively generated MHD waves with a one-dimensional plasma inhomogeneity, trans- 
verse to the magnetic field, is considered in the three-dimensional regime. Because of the transverse inhomogeneity, 
MHD fluctuations, even if they do not include initially any density perturbation, evolve toward states where the 
compressible components tend to become predominant. The propagating MHD pulse asymptotically reaches a 
quasi-steady state with the final levels of density perturbation weakly depending on the degree of non-planeness 
of the pulse in the homogeneous transverse direction and somewhat stronger depending on plasma 13. Our study 
demonstrates the necessity of incorporation of compressible and 3D effects in theory of Alfven wave phase mix- 
ing. However, as far as the dynamics of weakly non-plane Alfven waves is concerned it can still be qualitatively 
understood in terms of the previous 2.5D models. 
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1. Introduction 

Problems of heating of the open corona of the Sun and 
acceleration of the solar wind are closely related with the 
interaction of MHD waves with plasma inhomogeneities 
and, in particular, with coupling of compressible and in- 
compressible components of the waves. Investigation of 
coupling mechanisms between the waves is of prime im- 
portance because, on one hand, the Alfven waves are usual 
candidates of energy transport from the lower layers of 
the solar atmosphere to the corona, and, on the other, 
compressible perturbations are subject to much more ef- 
ficient dissipation than the incompressible Alfven waves 
(see, e.g., Narain & Ulmschneider (1996)). Also, the com- 



the dissipation is proportional to the wave number 
squared, phase mixing leads to enhanced dissipation of 
the Alfven wave directly. In the compressible plasma, 
as demonstrated by |Malara, Petkaki fc Veltri (1996) 



Nakariakov et al. (1997); 



Nakariakov et al. (1998) 



Botha et al. (2000); Tsiklauri, Arber fc Nakariakov (2001) 



pressive waves, as opposed to the Alfven waves, can trans- 
port energy across the magnetic field, due to the fact that 
their propagation in space is not constrained by the field 
lines. 

The origina l idea of the phase mixing o f incompressible 
Alfven waves ( Heyvaerts fe Priest (1983) ) was based on 
the following argument: when plasma has a density 
gradient perpendicular to the magnetic field, local Alfven 
speed is a function of the transverse coordinate. Thus, 
when an Alfven wave propagates along the field its per- 
turbations on the adjacent field lines become out of phase. 
This stretching of Alfven wave front creates progressively 
smaller spatial scales across the field. In turn, because 
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Tsiklauri, Nakariakov & Arber (2002), phase mix- 
ing of linearly polarized plane Alfven waves leads 
to the enhanced nonlinear generation of fast mag- 
netoacoustic waves. However, it was established, 
in 2.5D geometry, for the harmonic Alfven wave 
Botha et al. (2000) , and for a wide spectru m Alfven pulse 
by Tsiklauri, Arber fc Nakariakov (2001) that compres- 
sive perturbations, which are initially absent from the 
system, do not grow to a substantial fraction of the initial 
Alfven wave due the destructive wave interference effect. 

In this work we consider fully three dimensional 
geometry and study interaction of linear (with the 
non-linear effects totally ignored) MHD waves with a 
one-dimensional inhomogeneity of the plasma, taking into 
account compressibility of the plasma and the localization 
of the MHD pulse in the direction perpendicular to both 
the magnetic field and the inhomogeneity gradient. In 
particular, we shall study how the phenomenon of phase 
mixing is affected by these factors. A similar problem 
was studied by De Groof fc Goossens (2000) in the 
context of resonant absorption of MHD waves in coronal 
loops as a heating mechanism. However, our treatment 
is more general as we perform a direct numerical 3D 
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simulation without resorting to Fourier transform in time 
and y direction (and, consequently, our study can easily 
be generalized to the case of 2D and 3D structuring) 
allowing us to consider an initial value problem. Besides, 



Dc Groof & Goossens (2000) used a harmonic form of the 
initial perturbation as opposed to our spatially localized, 



Gaussian, one. Recently, Hood, Brooks & Wright (2002) 



showed that the wide spectrum regime of phase mixing 
can be quite different from the harmonic one, as in fact, 
phase mixing of localized Alfven pulses results in a slower, 
algebraic, damping as opposed to the standard exponen- 
tial damping of harmonic Alfven waves, which suggest 
that localized Alfvenic perturbations will transport energy 
higher into the corona than harmonic ones. Such Alfvenic 
perturbations could be generated e.g. by transient events 
such as solar flares, coronal mass ejections, etc. While 
most of the phase-mixing studies concentrated on har 
monic perturbations (e.g. |Hcyvaerts fc Priest (1983 



electron beam generation) of the Alfven waves in the 
auroral density inhomogeneities (cavities). 

In this work we study the propagation part of the prob- 
lem, i.e. we consider an ideal plasma limit. The work is in 
progress to include finite plasma resistivity in order to in- 
vestigate quantitatively the dependence of the decay of 
Alfvenic part of the MHD pulse upon the coupling to the 
existing compressive waves. 

The paper is organized as follows: in section 2 we for- 
mulate our model. In section 3 we present the results of 
numerical simulation, while we close in section 4 with the 
discussion of main results. 

2. The model 

In our model we use equations ideal MHD 



Hood, Ireland fc Priest (1997) 



Hood et al. (1997); 



Nakariakov et al. (1998) 



Ruderman, Nakariakov & Roberts (1998) 



DeMoortel et al. (1999); DeMoortel, Hood fc Arber (2000) 



Grappin, Leorat fe Buttighoffer (2000) 



Botha et al. (2000)), only few considered spa- 
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p— + p(V ■ V)V = -Vp B x curLB, 

dt Ait 

ri B 
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% + V ■ Vp + 7PV • V = 0, 
ot 



tially localized ones 


(Nakariakov et al. (1997) 


Tsiklauri, Arber & Nakariakov (2001) 




Tsiklauri, Nakariakov & Arber (2002) 


Hood, Brooks & Wright (2002)) 







An additional motivation to this study is connected 
with the growing interest to the problem of interaction of 
MHD waves with 2D and 3D plasma structures and irreg- 
ularities, as the coronal and wind plasmas are observed to 
be structured in all three dimensions. 2D and 3D structur- 
ing can dramatically affect properties of MHD waves. In 
particular, as it was show n in the incompressible regime by 
Similon fc Sudan (1989)| and confirmed in numerical ex- 
pcri mcnts performed by Petkaki, Malara fc Vcltri (1998) 



and [Malara, Petkaki fc Veltri (2000)| , the structuring dra- 
matically increases the efficiency of wave dissipation. 
However, in those studies, the compressible effects were 
not taken into account. From this point of view, our in- 
vestigation of the interaction of a 3D MHD pulse with a 
ID plasma inhomogeneity in the compressible regime pro- 
vides a necessary building element of the general theory 
of MHD wave interaction with plasma inhomogeneities. 

The present model is based upon the MHD de- 
scription of plasma and we focus on the problems 
that are of relevance to the heating of solar corona, 
acceleration and the dynamics of solar wind. While, 
similar studies exist that deal with, for instance, ter- 
restrial auroral applications. Namely, an alternative to 
MHD, more relevant, particle-in-cell simulations (using 
a guiding center implicit code for the electrons) have 



been performed ( 


Genot, Louarn & Le Queau (1999) 


Genot, Mottez & Louarn (2000) 


), which study the 


Genot, Louarn & Mottez (2001) 



(1) 

(2) 
(3) 



where B is the magnetic field, V is plasma velocity, p is 
plasma mass density, and p is plasma thermal pressure. In 
what follows we use 5/3 for the value of 7. 

We solve equations (|lj)-(||) in Cartesian coordinates 
(x, y, z). Note that as we solve a fully 3D problem we 
retain variation in the y-direction, i.e. (d/dy 7^ 0). The 
equilibrium state is taken to be an inhomogeneous plasma 
of density po{x) and a uniform magnetic field Bo in the 
z-direction. We consider a plasma configuration similar to 
the one investigated in [Malara, Petkaki fc Vcltri (1996) 



Nakariakov et al. (1997); 



Nakariakov et al. (1998) 



Botha et al. (2000)| ; [Tsiklauri, Arber fc Nakariakov (2001) 
i.e. the plasma has a one-dimensional inhomogeneity in 
the equilibrium density po(x) and temperature T (x). 
The unperturbed thermal pressure, po, is taken to be 
constant everywhere. 

Next, we do usual linearization of the Eqs. (Q)-((|) and 
write them in component form as following 



Po(x) 



dV x | dp B (dB x dB z 
dt dx 47r \ dz 

dB z 



dV y dp B Q 

Po{x) ^r + d-y + ^ 



dy 



dx 

dBy 

dz 



= 0, 



= 0, 



Po(x) 



dt 



dp 

dz 



0. 



dt 

dB z 

propagation and collisionless dissipation (via effective Qi 



dB x dV x 
dt dz 

^-Bo^ = 0, 



dz 



Bn 



dVx 
dx 



dVy 

dy 



= 0, 



(4) 

(5) 

(6) 
(7) 
(8) 
(9) 
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Fig. 1. Contour plots of V y and p'/po(x) (relative density perturbation) through three cross-sections at the same time 
snapshots. Here, (3 = 0.5, a y — 0.6, a z = 1.0. 



dp , ( 9V X dV y dV z 



dt ' V dx ' dy 



(10) i d u ~ + 4(*)) d 2 yy - c 2 A {x)d 2 zz ] V y (12) 

[(c 2 s (x) + c\{x)) d 2 xy ] V x - [c 2 s (x)d 2 yz ] V z = 0, 



It is useful to re-write Eqs. in a form of three _ g^gg j ^ _ [^(^2 j ^ 

coupled wave equations as following 



(13) 



[d 2 tt - (c 2 s (x) + c\{x)) d 2 xx - c\{x)d 2 zz \ v x (ii) 

- [(c 2 (x) + c\{x)) 8%] V y - [c 2 s {x)d 2 xz ] V z = 0, 



- [c 2 s [x)d* z ] V y = 0, 

where c A {x) = B / y/ A-k p {x) and c s (x) = y/"{Po/po{x) 
denote local Alfven and sound speeds respectively. 
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Fig. 2. A snapshot of relative perturbation of plasma den- 
sity by an MHD pulse interacting with a one-dimensional 
inhomogeneity at t — 15 through x — cross-section. 
Here, j3 — 0.5, a y — 0.6, a z = 1.0. 



We solve Eqs.(|4])-(|lCl|) numerically after re- writing 
them in a dimcnsionless form using following normal- 
ization: B 



x,y,z 



— BqB x 



(x,y,z ) = a*(x,y,z ), 

c A (x) = B Q /y/An PQ {x ) = B /y/Wp^/y/3 - 2 tanh(Aa;) = 
c \l \Z^> — 2 tanh(Ax), t = (a*/c A )i, V x , y , z = c* A V x , y , z . 
Note, that c s (x) — y^y/3/2c A (x) , where /3 stands for the 
ratio of thermal and magnetic pressures (3 = p / (B 2 / 8tt) . 
Here, A is a free parameter which controls the steepness 
of the density profile gradient. In our simulations we use 
A = 0.5. In what follows we omit bars on top of the phys- 
ical quantities. 



Fig. 3. The same as in Fig. (^), but through 
section. 



cross- 



(*))• 




3. Numerical Results 

In order to solve Eqs.(^)-(|l(i|) numerically we have writ- 
ten a new numerical code dtJ^dxlO, which uses a high 
order finite difference scheme. Namely, it evaluates 10- 
th order centered spatial derivatives, and advances so- 
lution in time using 4-th order Runge-Kutta algorithm. 
Therefore, dt4dxl0 is 0[(At/T) 4 ] accurate in time and 
0[(max(Ax, Ay, Az)/L) 10 } accurate in space, with T and 
L denoting run time and linear size of the simulation do- 
main. The use of such a high-order numerical approach 
was motivated by the very nature of the problem consid- 
ered, as development of phase mixing leads to the gener- 
ation of very steep profiles in the wave front. 

The simulation cube size is set by the limits —25.0 < 
x < 25.0, -25.0 < y < 25.0 and -25.0 < z < 25.0. 
Boundary conditions used in all our simulations are zero- 
gradient in all three spatial dimensions. 

We have performed calculation on various resolutions 
in attempt to achieve convergence of the results. The 
graphical results presented here are for the spatial res- 
olution 128 3 , which refers to number of grid points in x, y 
and z directions respectively. We have also performed cal- 
culation on the spatial resolution 256 3 and we found that 



Fig. 4. The same as in Fig. 
section. 



), but through y = cross- 



the results converge perfectly. This is understandable due 
to the high order of the scheme even for the 128 3 resolu- 
tion time-error, 0[(At/T) A ] sa 7 x 10~ 8 , and spatial -error, 
0[(max(Ax, Ay, Az)/L) w ] » 8 x 10~ 22 . Since a 256 3 res- 
olution run takes about 8 hours on Compaq ES40 with 
eight EV6 500-MHz processors, while a 128 3 resolution 
takes only 1 h on Compaq ES40 with four EV6 500-MHz 
processors producing the same results, we opted for the 
latter, less CPU-consuming, alternative. 

In the numerical simulations the MHD perturbation 
is initially a plane (with respect to ^-coordinate) pulse, 
which has a Gaussian structure in y and z-coordinates 

V y (x, y,z,t = Q) = c A (x) exp (-{a y y) 2 - (a z z) 2 ) . (14) 

Here, a y and a z are free parameters which control the 
strength of gradients in y and z direction of the initial per- 
turbation. As the problem considered is linear, the wave 
amplitude can be taken to be normalized to unity. 

In the geometry considered, when the initial perturba- 
tion depends on y it cannot be regarded as pure Alfvenic 
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Fig. 5. Maximum of the absolute values over the whole 3D simulation box of all physical quantities as a function of 
time. Here, f3 = 0.5 and a z = 1.0. Thick solid lines correspond to the resolution 256 3 , while thin solid lines to that of 
128 3 , both for a y = 0.6. Note that on all seven panels these two curves practically overlap, which serves as a proof of 
convergence of our simulation results. Dashed lines represent a y — 0.4, while dash-dotted lines show a y = 0. 



one. In fact, in such a pulse, all three waves - Alfven, 
fast and slow magnetosonic waves - are inter-coupled so 
that there is no use of their separation per se. The pulse 
is set to be initially plane in the x-direction. This allows 
us to emphasize the effect of the inhomogeneity on the 
pulse evolution. From the point of view of applications, 
this simply means that the initial characteristic size of 
the pulse in that direction is greater than the scale of the 
inhomogeneity. 

The particular choice of the initial condition Eq. ( |l4"| ) 
is motivated by the argument that when we set B y ini- 
tially to zero we automatically guarantee fulfillment of 
divi? = when a y ^ 0, and as the system evolves it 
adjusts itself which mode to excite (depending whether 
a y is zero or not). For instance, if we choose our ini- 
tial conditions as V y given by Eq. (|l4|), B y (x 1 y,z,t = 
0) = exp (— (a y y) 2 — (a^) 2 ) , with a y — and the 
rest of physical quantities set to zero, we would excite 
a pure Alfvenic pulse traveling in the negative direction 
along z-axis. However, as long as a y ^ we have to set 
B y (x, y, z, t = 0) = in order to fulfill divB = 0. Besides, 
our choice of the localized "kinematic" (the velocity is ini- 
tially perturbed only, while the magnetic field is constant) 



perturbation Eq. ([l4]) is well motivated by the fact that 
such perturbations can arise, e.g., during coronal mass 
ejections, solar flares, or other violent events (see, e.g., 
Roussev et al. 2001), etc. 

In Figure [l] we present time evolution of the transverse 
component of the plasma velocity, V y and the density per- 
turbation p'/pq(x). The choice of these physical quanti- 
ties is motivated by the following reasoning: the trans- 
verse velocity component V y is subject to Alfven wave 
phase mixing and the density perturbation demonstrates 
the effect of the plasma compressibility. The left column 
of Fig. [l| shows three snapshots of V y and two snapshots 
(p'/po(x)) through the y = cross-section, the middle 
column shows three snapshots of V y and two snapshots 
(p' / po(x)) through the x = cross-section, and the right 
column shows three snapshots of V y and two snapshots 
(p' /po(x)) through the z — cross-section. The snapshots 
of V y and p' / po(x) at times t = 7.5 and 15 are placed such 
that it is easier to spot the differences. The snapshot of 
p' / po(x) at t = is not present in the figure as it is iden- 
tically zero and it is generated later by the y-gradients in 
the function V y . This compressible perturbation is a sig- 
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Fig. 7. All physical quantities perturbed by the MHD pulse as a function of time at a given point in the phase mixing 
region. Here, j3 = 0.5, a y — 0.6, a z — 1.0. Left column: at the point (x=0,y=-8,z=0), mid column: at the point 
(x=0,y=-8,z=4), right column at the point (x=0,y=-8,z=8) 



nificant component of the pulse. What we gather from this 
graph is as following: 

— y = cross-section: the initial MHD pulse, which 
is plane with respect to x-coordinate, is split in two 
D'Alambert's solutions, with half of the amplitudes, 
traveling in two opposite directions along the magnetic 



field. Because the unperturbed density is inhomoge- 
neous across the magnetic field (in the rc-direction) the 
local Alfven speed depends upon x. Thus, the Alfvenic 
pulse is phase-mixed and its initially plane front is 
continuously distorted creating transverse gradients. 
Note, that the right wing of the pulse is asymmetric to 
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Fig. 8. All physical quantities as a function of time at a given point far from the phase mixing region. Here, (3 — 0.5, 
a y = 0.6, a z = 1.0. Left column: at the point (x=16,y=-8,z=0), mid column: at the point (x=16,y=-8,z=4), right 
column at the point (x=16,y=-8,z=8) 



the left one, this is due to the fact that its amplitude 
drops from 1.0 in the x > domain to l/y/E in the 
x < 0. The generated density perturbation evolves in 
the similar manner, however there is a notable differ- 
ence - as plasma [3 = 0.5 in this case, we observe that 
the density perturbation travels at a slower speed (slow 



magnetosonic wave speed) than the Alfvenic one. If we 
look at the right wing of V y (x,0, z) at t = 15, we see 
that it traveled to z=15, which is consistent with the 
Alfven speed equals 1 (recall that for x > ca(x) = 1). 
The density perturbation as it travels along the field 
should have a velocity of a slow magnetosonic wave, 
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Fig. 6. Normalized kinetic, magnetic, internal and total 
energies (all of them integrated over the volume) as a func- 
tion of time. Here, (3 = 0.5 and a z = 1.0. Solid lines cor- 
respond to a y = 0.6, while dashed and dash-dotted lines 
show cases of a y = 0.4 and a y = respectively. 



which for the field aligned propagation coincides with 
the local speed of sound. For the parameters consid- 
ered c s (x) — y/^P/2cA(x) = 0.65ca{x). Therefore, the 
right wing of the density perturbation at t = 15 should 
have traveled to a position z w 9.7. That is what we 
actually see from the bottom panel of the left column 
in Fig§ 

— x = cross-section: the initial Alfvenic pulse which 
has an asymmetry as a y ^ a z Gaussian bell shape is 
split in two smaller amplitude bells traveling in two 
opposite directions along the magnetic field plus some 
semi-elliptic shaped wake expanding outwards in all di- 
rections. Note, that because the unperturbed density 
is homogeneous across the magnetic field in y-direction 
we do not see any distortions on the wave front with 
respect to the y — line as we saw in the previous 
case. The generated density perturbation evolves in the 
similar manner, however a notable difference (can be 
better seen from Fig. ||) is that the solution consists of 
positive and negative spikes plus semi-elliptic shaped 
wakes that are symmetric with respect to p' / po(x) = 
plane. This can be understood readily as due to the 
conservation of total mass, integral of density pertur- 
bation over the volume should, indeed, be zero. 

— z = cross-section: the MHD pulse evolution, viewed 
through this cross-section, is similar to the case of 
y = cross-section (see above). Again we see two 
phase-mixed D'Alambert's solutions propagating in 
the opposite directions. The notable difference is that 
through z — cross-section the pulse looks wider as 
a y = 0.6 while a z = 1.0. The generated density pertur- 
bation evolves in the similar manner, as the Alfvenic 
one. However, in differ from y = cross-section, we 
see no difference in traveling speeds. This can be ex- 
plained as following: the pulse was not excited as 



a "pure Alfven wave". In fact, should the perturba- 
tion V y be purely Alfvenic there would have been no 
propagation along y axis at all (recall that the mag- 
netic field is directed along the z-axis). Thus, what 
we see through the z = cross-section both for V y 
and p' I po(x) is the compressive part of the perturba- 
tion traveling at the fast magnetosonic wave speed. 
For the parameters considered, the fast magnetosonic 
wave speed is a/1 + r y(3/2cA{x) — 1.2ca{x). Therefore, 
the right wing of both V y and density perturbation at 
t = 15 should have traveled to a position y w 17.9. 
That is what we actually observe from the two bottom 
panels of the right column in Fig.0. 

Also, the figure shows the obvious anisotropy in the pulse 
propagation, connected with the presence of two chosen 
directions: the direction of the straight magnetic field and 
the direction of the density gradient. 

Perturbation of B y is initially absent from the system, 
but as it is a potential component of the kinetic coun- 
terpart, Vy, of the wave, it is soon generated and further 
evolves similarly to V y in a form of two negative and posi- 
tive phase-mixed pulses traveling into two opposite direc- 
tions. Also, the development of the pulse is accompanied 
by generation of compressible components. The structure 
of the density perturbations in the pulse at the time t = 15 
is shown in Figures. ^, || and || Note, that in all the figures, 
pre-factor e stands to underscore the fact that the linear 
problem, which we study, does not depend on the initial 
amplitude. Thus, we use A = 1, while we have to bear 
in mind that linear approximation itself is valid for small 
amplitudes that is why we use small, arbitrary pre-factor 
e in our notations. 

A fairly good quantity describing the relation between 
different components of the pulse is the maximum of ab- 
solute value over the whole simulation domain. We plot 
this quantity for all physical variables in Figure ||. In this 
figure, thick solid lines correspond to the resolution 256 3 , 
while thin solid lines to that of 128 3 , both for a y = 0.6. It 
is remarkable that on all panels these two curves practi- 
cally do overlap, which serves as a proof of convergence of 
our simulation results. The last figure in the bottom row 
presents the maximum of absolute value over the whole 
simulation domain of div B , and we indeed observe almost 
perfect fulfillment of the fundamental law, div£? = 0, 
which comes as a bonus of high-order, centered, finite dif- 
ference numerical scheme. As expected, Fig. |^ shows that 
when a y = 0, the pulse is perfectly Alfvenic (perturb- 
ing V y and B y only). This conclusion can be deducted 
either from analyzing Eqs. (|ll|)-(|l3|) or resorting to a clas- 
sic mechanical analogy of coupled pendulums. In effect, 
Eqs. (pr|)-(p^|) also describe three inter-coupled mathe- 
matical pendulums which are located in xOz plane. Thus, 
as long as we do not perturb these pendulums such that 
a y ^ 0, they will always oscillate in the xOz plane. Thus, 
what we see in Fig. || when a y = (dash-dotted lines) 
is that initial kinematic Alfven perturbation (V y ) is split 
in half- amplitude D'Alambert's solutions (both for V y and 
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By) and no other physical quantity is generated. However, 
when we switch on the coupling, a y 7^ 0, compressible per- 
turbations are generated. 

In order to investigate the effect of phase mixing on the 
created quasi-steady MHD state, we produced time series 
of all physical quantities in several points of the simulation 
3D box, Figures ^ and |^. Namely, at the points (x=0,y=- 
8,z=0, 4, 8), which are located in the phase mixing (spatial 
inhomogcneity) domain, and at (x=16,y=-8,z=0, 4, 8), 
which are far away from it. The three panels (from the 
left to the right), in effect, trace the dynamics in the z- 
direction, i.e. along the regular magnetic field lines. These 
types of data would be obtained by, for example, three 
satellites which are located in three different regions of 
solar wind. 

There are two noteworthy features that can be gath- 
ered from these plots: First, we see that in the phase mix- 
ing region, V x , which is associated with the fast magne- 
tosonic component of the pulse, attains about 40 times 
larger values than in the region that is far from the phase 
mixing region. That is a sensible result, since it is known 
that phase mixing efficiently generates oblique fast magne- 
tosonic waves. Second, the total magnetic field perturba- 
tion \B — .Bo I is positively correlated with the pressure (as 
well as density, which is proportional to the gas pressure) 
perturbation, which indeed is what was expected from the 
fast magnetosonic wave. 

Also, Figs demonstrate that, interacting with the 
plasma inhomogeneity, the longitudinally (parallel to the 
magnetic field) propagating part of the initial perturba- 
tion develops to an almost non-evolving (without change 
in amplitude) propagating state. In this quasi-steady 
MHD state all components of the pulse tend to propa- 
gate without change in amplitude, even though the effect 
of phase mixing is in action all the time (^-component of 
the inverse characteristic spatial scale tends to infinity). 

In Figure ^ we investigate the energetics of our numer- 
ical simulation. In fact, we observe nearly perfect (±0.2% 
error) conservation of total energy by dt^dxlO numerical 
code. The major conclusion which can be drawn from this 
graph is that when a y — there is no internal (compres- 
sive) energy generation, while with the increase of a y its 
final, asymptotic levels increase progressively. 

The parametric space of the problem is studied in 
Figure ^ by plotting the final levels of relative density per- 
turbation, approached by the MHD pulse, as a function of 
the initial localization of the pulse in the y-direction a y , 
and plasma (3. We gather from Fig. |^ that the achieved 
levels of density fluctuations depend weakly on a y , while 
there is somewhat stronger dependence on (3. 

4. Conclusions 

In this paper we present results of numerical modeling 
of interaction of an initially localized MHD pulse with a 
transverse inhomogeneity of the plasma. The pulse is non- 
plane in all three spatial directions and is subject to the 
effect of phase mixing. An impulsively generated perturba- 




0.2 0.4 0.6 0.0 1.0 1.2 1.4 1.6 



Fig. 9. Top panel: the final levels of relative density per- 
turbations as a function of a y (for (3 = 0.5 and a z = 1.0). 
Bottom panel: the same but as a function of (3 (for 
a y = 0.6 and a z = 1.0). 

tion of the plasma velocity develops to an anisotropically 
propagating pulse, which has a significant compressible 
component. More specifically: 

— The non-uniformness of the pulse in the homogeneous 
transverse (or, in other words, in the third, perpendic- 
ular to both the magnetic field and the inhomogeneity 
gradient) direction (non-zero a y ) leads to appearance 
of compressive perturbations in the pulse. 

— In the presence of plasma inhomogeneity, a non- 
uniform MHD pulse is essentially compressible. More 
specifically, when the pulse is initially plane in the ho- 
mogeneous transverse direction, a y = 0, there is no 
internal (compressive) energy generation, while with 
the increase of a y , the levels of the compressible en- 
ergy increase progressively. 

— A propagating MHD pulse asymptotically reaches a 
quasi-steady state. The final levels of density perturba- 
tion, which can be considered as a measure of the com- 
pressibility in the pulse, depend weakly on a y , while 
there is somewhat stronger dependence on (3. 

— A smooth ID transverse inhomogeneity of the plasma 
supports propagation of compressible MHD pulses. 
This mechanism of wave guiding is different from the 
well-understood phenomenon of refraction of fast mag- 
netosonic waves, as the profile of the Alfven speed in 
the inhomogeneity does not have a minimum. 

— Weakly non-plane MHD waves are subject to phase 
mixing and can qualitatively be considered in terms 
of the 2.5D theory of Alfven wave phase mixing. 
Consequently, the 2.5D theory of Alfven wave phase 
mixing remains relevant to the 3D case too. However, 
quantitative theories of the interaction of MHD waves 
with plasma inhomogeneities should include the com- 
pressibility of the plasma as a necessary ingredient. 

Our main conclusion is that the effects of three- 
dimensionality, compressibility and inhomogeneity should 
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be all together taken into account in the wave-based the- 
ories of coronal heating and solar wind acceleration, and 
as well as in the theories of MHD turbulence. 

Acknowledgements. The authors are grateful to Tony Arber 
for valuable advise when dt^dxlO code was written and for 
a number of valuable discussions. DT acknowledges financial 
support from PPARC. Numerical calculations of this work were 
performed using the PPARC funded Compaq MHD Cluster at 
St Andrews and Astro-Sun cluster at Warwick. 

References 

Botha, G. J. J., Arber, T. D., Nakariakov, V. M., & Keenan, 

F.P 2000, A&A, 363, 1186 
De Groof, A., & Goossens, M. 2000, A&A, 356, 724 
DeMoortel, I., Hood, A. W., Ireland, J., & Arber, T. D. 1999, 

A&A, 346, 641 

DeMoortel, I., Hood, A. W., & Arber, T. D. 2000, A&A, 354, 
334 

Genot, V., Louarn, P., & Mottez, F. 2001, J. Geophys. Res., 
106, 29633 

Genot, V., Mottez, F., & Louarn, P. 2000, J. Geophys. Res., 
105, 27611 

Genot, V., Louarn, P., & Le Queau, D. 1999, J. Geophys. Res., 
104, 22649 

Grappin, R., Leorat, J., & Buttighoffer, A. 2000, A&A, 362, 
342 

Heyvaerts, J., & Priest, E. R. 1983, A&A, 117, 220 

Hood, A. W., Brooks, S. J., & Wright, A. N. 2002, Proc. Roy. 

Soc. A (in press) 
Hood, A. W., Ireland, J., & Priest, E.R. 1997, A&A, 318, 957 
Hood, A. W., Gonzalez-Delgado, D., & Ireland, J. 1997, A&A, 

324, 11 

Malara, F., Petkaki, P., & Veltri, P. 2000, ApJ, 533, 523 
Malara, F., Primavera, L., & Veltri, P. 1996, ApJ, 459, 347 
Nakariakov, V. M., Roberts, B., & Murawski, K. 1997, Solar 
Phys., 175, 93 

Nakariakov, V. M., Roberts, B., & Murawski, K. 1998, A&A, 
332, 795 

Narain, U., & Ulmschneider, P. 1996, Space Sci. Rev., 75, 453 
Petkaki, P., Malara, F., & Veltri, P. 1998, ApJ, 500, 483 
Ruderman, M. S., Nakariakov, V. M., & Roberts, B. 1998, 

A&A, 338, 1118 
Roussev I., Galsgaard K., Erdelyi R., Doyle J.G. 2001, A&A 

370, 298 

Similon, P. L., & Sudan, R. N. 1989, ApJ, 336, 442 
Tsiklauri, D., Arber, T. D., & Nakariakov, V. M. 2001, A&A, 
379, 1098 

Tsiklauri, D., Nakariakov, V. M., & Arber, T. D., 2002, A&A 
(submitted) 



